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Abstract 

We report a measurement of the ti production cross section using the CDF II detector at the 
Fermilab Tevatron. The data consist of events with an energetic electron or muon, missing trans- 
verse energy, and three or more hadronic jets, at least one of which is identified as a 6-quark jet 
by reconstructing a secondary vertex. The background fraction is determined from a fit of the 
transverse energy of the leading jet. Using 162±10 pb -1 of data, the total cross section is found 
to be 6.0±1.6(stat.)±1.2(syst.) pb, which is consistent with the Standard Model prediction. 

PACS numbers: 14.65.Ha, 13.85.Ni, 13.85.Qk 
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I. INTRODUCTION 



The top quark is the most massive of nature's building blocks yet discovered. Because new 
physics associated with electroweak symmetry breaking will likely couple to an elementary 
particle in proportion to its mass, it is important to measure the top quark couplings as 
accurately as possible. In the strong interaction sector, the couplings are reflected in the 
it production cross section in hadron collisions. Previous measurements were made in pp 
collisions at a center-of-mass energy of 1.8 TeV We recently reported a result Q|, using 
data taken at 1.96 TeV with the CDF II detector at the Tevatron collider, using the double 
leptonic decay mode of the top quark. Here we report a measurement of the it production 
cross section using a different decay mode and a new method. 

In order to measure the cross section, one first has to obtain a sample rich in top quarks 
and then determine the amount of background in the sample. We select events consistent 
with the decay chain it — ► WbWb — > lubqq'b, where the charged lepton I is either an 
electron or muon. We start with events containing an energetic electron or muon, significant 
transverse momentum imbalance indicative of a non-interacting neutrino, and at least three 
hadronic jets. To enrich the sample in top quarks, we require that at least one jet contain 
a secondary vertex consistent with the decay of a B hadron. 

Measurements in the past have relied on the ability of theoretical calculations to deter- 
mine, for the background, the fraction of events that contain 6-quark jets. In this paper, 
we instead measure the background fraction directly in the signal data sample. The trans- 
verse energy of the highest Et jet j^| or the second highest E? jet is a good discriminator 
between signal and background. Typically in a tt event, these jets are the primary decay 
products (6-jets) of the very heavy top quarks and thus have a hard Et spectrum. For most 
of the background sources, however, they are produced as QCD radiation, resulting in a 
much softer bremsstrahlung-like Et distribution. We use the leading jet Et spectrum for 
the primary measurement of the signal fraction. The second leading jet distribution is used 
to check the result. 

In order to obtain the background spectrum, we need data that are kinematically similar 
to our final sample, but which do not have significant tt contamination. We show that the 
leading jet Et spectra for the background processes are similar whether or not the events 
contain 6-quark jets. Then the non-heavy flavor spectrum becomes the background template 
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for measuring the tt fraction in the signal 




We use the HERWIG |J and PYTHIA 



Monte Carlo calculations followed by a simu- 



lation of the CDF II detector to obtain the tt signal behavior. The soft spectrum of the 
parton showers in these Monte Carlo models is not relevant for the signal shape. To test 
that our method is plausible, we study the background shape using the ALPGEN+HERWIG 
Monte Carlo ALPGEN provides a harder and more realistic jet Et spectrum. For the 
study of 6-jet identification, the PYTHIA calculation is used. 

II. DETECTOR 

The CDF II detector |?| is an azimuthally and forward-backward symmetric apparatus 
designed to study pp collisions at the Fermilab Tevatron. It consists of a magnetic spectrom- 
eter surrounded by calorimeters and muon chambers. The charged particle tracking system 
is immersed in a 1.4 T magnetic field parallel to the p and p beams. A 700,000-channel 
silicon microstrip detector (SVX+ISL) provides tracking over the radial range from 1.5 to 
28 cm. A 3.1 m long open-cell drift chamber, the Central Outer Tracker (COT), covers the 
radial range from 40 to 137 cm. The COT provides up to 96 measurements of the track 
position with alternating axial and 2° stereo superlayers of 12 wires each. The fiducial region 
of the silicon detector extends to ~ 2, while the COT provides coverage for < 1. 

Segmented electromagnetic and hadronic calorimeters surround the tracking system and 
measure the energy of interacting particles. The electromagnetic and hadronic calorime- 
ters are lead-scintillator and iron-scintillator sampling devices, respectively, covering the 
pseudorapidity range |?7| < 3.6. The electromagnetic calorimeters are instrumented with 
proportional and scintillating strip detectors that measure the transverse profile of elec- 
tromagnetic shower candidates at a depth corresponding to the shower maximum. Drift 
chambers located outside the central hadron calorimeters and behind a 60 cm iron shield 
detect muons with \rj\ < 0.6. Additional drift chambers and scintillation counters detect 
muons in the region 0.6< \r]\ < 1.0. Gas Cherenkov counters measure the average number 
of inelastic pp collisions and thereby determine the luminosity with the coverage 3.7< < 



The results reported here are based on data taken in Fermilab Collider Run II between 
March, 2002 and September, 2003. The integrated luminosity is 162 pb _1 for events selected 



4.7. 



9 



with an electron or central muon. For muon events with \r]\ between 0.6 and 1.0, the 
integrated luminosity is 150 pb _1 . 

III. EVENT SELECTION 

A. Lepton trigger 

CDF employs a three level trigger system, the first two consisting of special purpose 
hardware and the third a farm of computers. For the electron top sample, the level-1 
trigger requires a track of Pt > 8 GeV/c matched to an electromagnetic calorimeter cell 
containing > 8 GeV with a small amount of energy in the hadronic cell behind it. 
Calorimeter energy clustering is done at level-2, and the > 8 GeV/c track must be matched 
to an electromagnetic cluster with above 16 GeV. At level-3, a reconstructed electron 
candidate with E? > 18 GeV is required. The Pt > 18 GeV/c level-3 muon triggers come 
directly from two level-1 triggers: a track with Pt > 4 GeV/c is matched to a stub in 
the central muon chambers; or a track with Pt > 8 GeV/c is matched to a stub in the 
0.6 < \rj\ < 1.0 muon chambers. 

B. W + jets selection 

After full event reconstruction, we require lepton candidates to pass identification criteria 
and to be isolated from other energy deposits in the calorimeter. The event selection criteria 
are the same as those in Ref. |8j, where they are described in detail. Electron candidates 
must have a well-measured track pointing at a cluster of energy in the calorimeter with Et 
> 20 GeV. The lateral and transverse shower size in the calorimeters as well as the trans- 
verse profile in the shower-maximum detectors must be consistent with an electromagnetic 
cascade. Muon candidates with P T > 20 GeV/c must pass through calorimeter cells whose 
energy deposition is consistent with the ionization of a muon, and the reconstructed posi- 
tion of the track segment in the muon chambers is required to be consistent with multiple 
Coulomb scattering of the extrapolated track from the COT. Lepton candidates must also 
be isolated. Isolation (/) is defined as the ratio between calorimeter energy in a cone of 
radius 0.4 in the rj-cp plane around the lepton, but excluding the lepton, divided by the 
lepton energy. We require / < 0.1. In addition, all candidate events must have $t > 20 
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GeV. The $t is corrected for both muon momentum and the position of the pp collision 
point. Jets are found using a fixed-cone algorithm with a cone radius of 0.4 in rj-(p space. To 
obtain the correct jet energy, this analysis applies three corrections after jet clustering. We 
correct for detector response variations in 77, detector stability, and a correction for multiple 
interactions in an event. For this analysis, jets are counted if they have > 15 GeV and 
\rj\ < 2.0 after all the corrections are applied. We select events with three or more jets to 
retain high acceptance for ti events, allowing one jet to fail our Ej> or rj requirement. 

C. b-jet identification 

In this sample, the major background to ti is the electroweak production of a W boson 
with hadron jets produced by QCD. These W + jets events usually contain only light quark 
and gluon jets, whereas signal events always contain two b- quark jets. Thus identification 
of 6-jets (6-tag) provides a significant increase in the signal-to-background ratio. 

We identify 6-quark jets through the metastable B hadrons in the jet fragmentation. 
Their ~ 1 ps lifetime translates into a secondary vertex a few millimeters from the primary 
interaction. We use the excellent position resolution of the SVX+ISL to find these secondary 
vertices. The algorithm [9] proceeds as follows: (1) select at least two good tracks in a jet 
with both COT and SVX+ISL information, (2) search for a high quality secondary vertex 
using the selected tracks, (3) measure the distance in the transverse plane {Lr) between the 
primary and the secondary vertices, and (4) accept the secondary vertex if L T /a(L T ) > 3, 
where a(L T ) is the L T resolution. 

Based on simulation of the 6-tagging algorithm, we determined that requiring at least 
one of the jets in an event to be tagged as a b-jet is expected to retain 53% of the top quark 
events while removing more than 95% of the background events |9(. The measured cross 
section depends on the value of the 6-tagging efficiency used to extract it. The difference 
between the efficiency in the simulation and that in data has been measured with a 6-enriched 
sample of dijet events in which an electron is found in one jet (e-jet) and a secondary vertex 
is found in the other jet. From the fraction of e-jets that have an observed secondary vertex, 
we find that the simulation has a 6-tag efficiency higher than the data by 21%, a factor 
that is independent of jet E? |9j. This difference is corrected for in section VI A and the 
uncertainty is discussed in section VI B. 
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IV. BACKGROUND DETERMINATION 



There are several background sources in the 6-tagged I +$t + >3 jets sample. We 
give their approximate contributions to the total background, obtained from both data and 
theoretical calculations, although the precise composition is not needed for this analysis 
because the spectra are very similar to each other. The sources are the production of a W 
boson accompanied by the QCD production of heavy flavor quarks (W+HF: ~ 40% of the 
total number of the background events in W+ > 3 jets), mis-identification in 6-jet tagging 
mainly due to track and vertex resolution (mistag: ~ 30%), a fake W boson associated with 
a real or fake lepton (jaoia-W: ~ 20%), dibosonproduction (WW/WZ/ZZ) and electroweak 
tb production (diboson + single top: ~ 10%) 9|. 

A. Method 

We find that the leading jet E? is the best discriminant between the signal and back- 
ground among single jet Et variables after considering both statistical and systematic effects, 
including the difference between Leading-Order (LO) and a Next-to-Leading-Order (NLO) 
simulation [l^|. Consequently we use the shape of the leading jet Et spectrum to determine 
the signal and background fractions. For the background, we extract from data the shape for 
all of the sources except for the small diboson and single top components. We assume that 
the VF+HF and mistag shapes are the same, so that we can extract that distribution from 
the background-rich sample of events in which no jet is 6-tagged. To avoid subtle kinematic 
differences based on jet rapidity, we still require that at least one jet is taggable, i.e., has at 
least two good tracks that pass through the SVX+ISL detector. The assumption that the 
leading jet Et in W + multijet events is independent of the flavor content of the jets is first 
studied in Monte Carlo simulations and then tested using events containing a W boson and 
either 1 or 2 jets (see section IV B). 

Figure ^ shows the shape comparison of the highest jet Et between W+ three light 
flavor jets (W+LF) and the various heavy flavor contributions as predicted by ALP- 
GEN+HERWIG Monte Carlo calculation followed by a simulation of the CDF II detector. 
In ly+HF cases, at least one jet is required to be 6-tagged, while for H^+LF, we require that 
at least one jet is taggable. We then correct the W^+LF shape for the slight Et dependence 
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Leading jet Ey (GeV) 

FIG. 1: The leading jet Et spectrum for W plus three light flavor jets (solid histogram in all 
plots) compared to Wbb plus one light quark (left), Wc plus two light quarks (center), and Wcc 
plus one light quark (right). These are ALPGEN+HERWIG Monte Carlo calculations followed by 
a simulation of the CDF II detector, (normalized by area) 

in the 6-tag efficiency, which is taken from simulation (Fig. |2J). 

The effect of the 6-tagging algorithm could be different for IV+HF, where a real b-jet is 
tagged, and VF+LF, where a fake tag is found. However the Et dependences of the tagging 
efficiency and mistag probability are similar. Figure El compares VF+HF events with a 6-tag 
to PF+LF after applying the mistag probability measured in jet data. The agreement is 
good. 

The other large background comes from events that do not contain a W boson. The lepton 
candidate is either misidentified or due to semi-leptonic decay in a b- or c-jet. Such leptons 
are typically not isolated in the calorimeter. Consequently the shape of the leading jet Et 
spectrum for non-VF events is determined from events that still have large missing transverse 
energy {fir > 20 GeV), but whose lepton is not isolated. We assume that the leading jet Et 
shape of this sample is the same as for the non-PF background events in the signal region 
since lepton isolation is not correlated with the Et of other jets in an event. The non-VF 
background distribution is added to the other backgrounds with a relative normalization 
taken from absolute estimates of the various background sources. However, since the non- 
W jet Et spectrum is very similar to those in the jy+HF and mistag backgrounds, the final 
result is insensitive to the non-F/ fraction. A systematic uncertainty is taken based on a 
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FIG. 2: The efficiency of 6-jet identification as a function of the leading jet Et derived from 
ALPGEN Wbb Monte Carlo in the W+ > 3-jet sample. The curve is a fourth-degree polynomial 
fit below 140 GeV. 
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FIG. 3: The leading jet Et spectrum for mistagged W plus three light flavor jets (solid histogram 
in all plots) compared to Wbb plus one light quark (left), Wc plus two light quarks (center), and 
Wcc plus one light quark (right) in VF+3 or more jets sample. For W^+HF sample, we require at 
least one tagged 6-jet. Note that the mistag prediction is realistic since it is obtained from jet data. 
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large variation in this background fraction. 

The spectra from diboson and single top production are estimated from ALPGEN+ 
HERWIG Monte Carlo calculations. The leading jet Et spectrum from single top production 
is added to the total background shape using the theoretical cross section, which is 6% of the 
total background, while the diboson spectra are neglected because these spectra are similar 
to the other dominant background sources, and their contribution is expected to be small, 




B. Test using W + 1-jet and W + 2-jet data 

The performance of the background modeling is tested using events containing a high 
Pt lepton, large $t, and either 1 or 2 jets (recall the signal sample has 3 or more jets). 
This sample contains all of the signal region background sources, but the ti contribution 
is small, only ~ 4%. A similar procedure is applied here as in the W+ > 3-jet sample. 
The dominant background shape is taken from events without a 6-tag, but with at least one 
taggable jet, and then a non-W contribution is added with a fraction (~ 15.2%) determined 
from absolute background estimates ^. A correction is made for the Et dependence of 
the 6-tagging efficiency. The resulting spectrum should agree well with that of the 6-tagged 
events if our method is valid. Figure 0] shows the results using 309 jets in the W+l and 2 jet 
samples. The agreement is good, with a Kolmogorov-Smirnov (KS) test probability of 18%. 
This value is not the KS probability itself but the p-value based on pseudo-experiments 
using the maximum difference of the accumulating distributions. 



C. Background shape in the W+ > 3-jet sample 

The Et spectra of the four main backgrounds in the signal sample, W+ >3 jets, are 
shown in Fig. |3J The dominant contribution is taken from the non-6-tagged events (H^+HF, 
mistag). These spectra have been corrected for the shape of the b-tag efficiency, which 
has been determined from simulation, and is shown in Fig. El as a function of the jet Et- 
The small ti contamination, ~ 6%, in the non-6-tagged sample is subtracted iteratively, as 
described below, so that the amount of this contamination is consistent with the final ti cross 
section. The non-W component is 21% of the total and is shown as the shaded portion. The 
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FIG. 4: The leading jet Et distribution for the 309 jets in the W+1 and 2 jet data sample. 
The open histogram is the total background prediction, with the dark shaded region the non-W 
component. The data points are the 6-tagged events with statistical error bars. 

small contribution from single-top production is then added, and the diboson components 
(WW/WZ/ZZ) are neglected as noted above. 

As indicated above, the shape change due to the application of the 6-tagging algorithm is 
applied to the background spectra. The efficiency of the algorithm is defined as the number 
of events that have at least one 6-tagged jet divided by the number having at least one 
taggable jet. This function, which drops at very low Et, is fit to a fourth-degree polynomial 
below 140 GeV and a flat line above 140 GeV. Possible variations in this shape are considered 
as a systematic uncertainty on the cross section measurement. Note that only the shape, 
not the absolute efficiency, is used here. 

The final background shape is shown in Fig. |H1 We fit this shape for use in the final 
unbinned log-likelihood fit using a Landau distribution plus a Gaussian function. The fitted 
parameters are summarized in Table HI and the fit result is shown in Fig. El The y 2 /d.o.f. 
for the agreement between the fit function and the data points is 11.4/10. 
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FIG. 5: The calculated dominant background shape (VF+HF, mistag, non-PF) in VF+> 3 jets 
sample. The shaded part is the non-VF portion. Error bars are statistical only. 
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FIG. 6: The calculated leading jet Et distribution for the background in the W+ > 3-jet sample. 
The fitted curve shows the Landau plus Gaussian function that is employed in the final unbinned 
likelihood fit. 
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TABLE I: Presented in this table are the fit parameters used to describe the background prob- 
ability density function shown in Fig. HO (L) and (G) refer to the Landau function and Gaussian 
parameters, respectively. The variable MPV represents the most probable value of Landau func- 
tion. The means and sigmas are expressed in GeV. 

parameters values 
height(G) 30.80±10.27 
mean(G) 66.71±7.4 
sigma(G) 32.77±5.6 
height(L) 512.8±79.2 
MPV(L) 37.0±3.6 
sigma(L) 7.67±1.5 

V. SIGNAL FRACTION 

We use a HERWIG Monte Carlo calculation followed by a full detector simulation to 
obtain the tt signal shape. Figure [7| shows the predicted leading jet E T distribution for 
tt events. It is significantly harder than the background spectrum, making it possible to 
separate the two contributions using a fit to the data. The spectrum in Fig. is fit to a 
Landau distribution plus two Gaussians (x 2 /d.o.f. = 19.9/28). The fitted parameters are 
shown in Table ITT1 

To fit the data to a sum of the signal and background templates, we use an unbinned 
likelihood fit with the following form. 

N 

C = ]JP(E Ti ;R) 

i=l 
N 



[ J [R -^signal {Exi) + (1 — -R) -Pback ground (En) 



i=l 



where the signal fraction R=jf. — l+ ^ b na ' k is the one free parameter in the fit, P S i gn ai(ETi) 

is the signal probability density as a function of E?, and ^background (-Sri) is that of the 
background. We tested the ability of this fit procedure to report correct values of the signal 
fraction and its uncertainty using a large number of Monte Carlo pseudo-experiments. Each 
pseudo-experiment used the number of signal and background events in our data sample 
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TABLE II: Fit parameters for the tt leading jet Ex distribution in Fig. [7|using two Gaussians and 
a Landau distribution. (L) and (G) refer to the Landau and Gaussian parameters, respectively. 
The variable MPV represents the most probable value of Landau function. The means and sigmas 
are expressed in GeV. 



parameters 


values 


height(Gl) 


200.8±26.0 


mean(Gl) 


60.0±5.0 


sigma(Gl) 


35.6±1.1 


height(G2) 


-109.9±8.3 


mean(G2) 


-76.1±6.5 


sigma(G2) 


192.9±14.3 


height (L) 


3913.5±192.8 


MPV(L) 


65.8±1.1 


sigma(L) 


13.9±0.6 
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and used a range of signal fractions (R) centered around the value found from our fit to the 
data. 

As mentioned above, there is a small tt contamination in the untagged data sample that 
is used to create the background template. The amount of tt that is subtracted when making 
the template is determined by an iterative process. Initially the fit is done without removing 
a tt component from the background template. The number of top events reported by the 
fit is used along with the 6-tagging efficiency to calculate the number of top events in the 
untagged sample. A tt subtraction in the background template is then made, and the data 
are refit. This tt contamination is determined to be small, ~ 6%, thus only one iteration is 
necessary. The final background template after the iteration is shown in Fig. |U1 

The result of the fit of the W+ > 3-jet data sample is shown in Fig. |SJ The histogram 
contains the 57 data events in which at least one jet has been tagged as a fe-jet. The solid 
curve is the best fit, with the individual components shown as dashed (tt) and dot-dashed 
(background) curves. The insert contains — ln(£/£ 

max ) 9jS 3j function of signal fraction. The 

signal fraction obtained is R =0.68^Qjg. 

Although we selected the leading jet Et as the fit variable a priori, we have studied other 
variables to check the robustness of the result. For the second leading jet and the sum 
of the first and second leading jet E^s, we find signal fractions of R =0.751q;}3 (Fig. EJ) and 
R =0.65^Q;}g (Fig. Unj), respectively. The agreement is good. 



VI. tt CROSS SECTION 

A. Acceptance and Efficiency 

The ti cross section is obtained from the formula, 

f+T\ ^obs -Rfit /-, \ 

where A Q b s is the number of candidate W+ > 3-jet events with at least one 6-tagged jet 
(57 events), i?g t is the signal fraction determined from the likelihood fit (0.68^g;}g), and 
A t i is the geometric acceptance for tt events in the CDF II detector [8]. Note that this 
acceptance includes the branching ratios. The parameter e t i is the detector efficiency for tt 
events jsj, which includes the trigger, event vertex position, event 6-tagging, and the lepton 
identification efficiencies. It also includes the effects due to photon conversion, cosmic ray, 
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FIG. 8: The fifty-seven candidate events (histogram) with the best fit curve (solid). The 
best fit composition, ti (dashed) and background (dot-dashed), is also shown. The insert shows 
-ln(£/£ max J 

function of the signal fraction. 

dilepton, and Z° boson removal. The quantity J Cdt is the integrated luminosity. The term 
A t f€ t t was determined from a PYTHIA Monte Carlo [5] calculation and detector simulation 
with a number of individual efficiency components determined from the data. The result 
for Ae is 4.02±0.03(stat.)±0.43(syst.)%. The electron (muon) channel contributes ~57% 
(~43%) of the total Ae. All calculations have been done using a top quark mass of 175 
GeV/c 2 Q. Multiplying Ae by the integrated luminosity gives the denominator for equation 
(1), 6.42±0.8 pb- 1 . 
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2nd Leading jet E T (GeV) 



FIG. 9: The fit result using the second leading jet. The fifty-seven candidate events (histogram) 
with the best fit curve (solid). The best fit composition, ti (dashed) and background (dot-dashed), 
is also shown. The insert shows — ln(£/£ 

max J &- function of the signal fraction. 

B. Systematic uncertainties 

There are a number of sources of systematic uncertainty as summarized in Table IIHI 
Template shape uncertainties affect the signal fraction determination, while other effects 
mostly impact the acceptance. Systematic uncertainties in the signal fraction are determined 
by a series of pseudo-experiments in which the generated pseudo-data are changed based on 
the systematics and then refit using the original templates. If the systematic uncertainty 
affects both the template shapes and the acceptance the uncertainty is taken to be 100% 
correlated. 

The largest uncertainty originates from the effect of the jet energy scale on the tt sim- 
ulation. This comes from a number of sources including modeling the relative calorimeter 
response as a function of rj, the absolute hadron energy scale, the underlying event contri- 
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FIG. 10: The fit result using the sum of the leading and the second leading jets. The fifty-seven 
candidate events (histogram) with the best fit curve (solid). The best fit composition, ti (dashed) 
and background (dot-dashed), is also shown. The insert shows — ln(£/£ max ) as a function of the 
signal fraction. 

bution, and jet fragmentation 8]. The largest contributions are due to the rj correction of 
the jet energy and the energy scale uncertainty. The mean energy of the leading jet from 
top quark decay is varied by ±6.1%, or about 5 GeV, and this effect contributes 15.3% to 
the final top cross section uncertainty. The jet energy scale uncertainty does not contribute 
to the background template shape systematic uncertainty largely because it is determined 
from the data. 

There are uncertainties in both the absolute value and Et dependence of the 6-tag ef- 
ficiencies, which are determined from fe-jet rich and generic-jet control samples jjj. The 
uncertainty of the absolute 6-tagging efficiency is dominated by several sources: statistics of 
the control data and Monte Carlo sample, composition uncertainties of the control sample 
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and the branching ratio of the b semi-leptonic decay. The ratio of the 6-tag efficiency be- 
tween the Monte Carlo and 6-jet rich control data sample is formed, and we vary the value 
within the uncertainties to determine the change in cross section. The dependence of the 
6-tag efficiency uncertainty is determined using the slope difference of two ratios. We form 
as a function of E? the 6-tag efficiency ratio from the generic-jet data and the Monte Carlo 
simulation (slope and uncertainty), and also from the 6-jet rich data sample and Monte 
Carlo simulation. A weighted average of the two slopes is used to determine the overall 
slope uncertainty, which is then applied to the Monte Carlo simulation to determine the top 
quark cross section uncertainty. 

The uncertainty in the background shape due to the statistics of the background sample 
is estimated with a series of pseudo-experiments in which the contents of each bin in Fig. El 
is varied independently according to Poisson statistics, with the resulting distribution refit 
to get a new background spectrum. The luminosity uncertainty comes predominantly from 
the uncertainty in the total inelastic cross section. Uncertainties in the lepton identification 
efficiency, which affect the acceptance, are determined from events using the unbiased tracks 
of Z — > 11 decays in events with multiple jets (slO]. There are several other efficiency uncer- 
tainties due to the trigger efficiency, the photon conversion veto efficiency, the cosmic ray 
veto efficiency, and track finding efficiency. These systematic uncertainties are summarized 
in Table ITTT1 

The parton distribution functions for (anti-)protons affect not only the shape of the ti 
signal, but also the acceptances. These uncertainties are estimated by varying a s and the 
parton distribution functions within the universal fit uncertainty. There are also uncer- 
tainties from the amount of gluon radiation in the Monte Carlo generators. The amount 
of initial state radiation is studied using high mass Drell-Yan dilepton data. The non-W 
contribution to the background shape has an uncertainty both in its relative amount (to 
the overall background) and its shape. The former is estimated by a ±100% variation in 
the amount of background measured from the non-isolated lepton sample. The shape un- 
certainty is measured by changing the nominal mixture of events containing non-isolated 
electron and muon candidates (~ 3.7 : 1 electrons to muons) in the data control sample to 
either 100% electrons or 100% muons. The shapes are used in the fit and the change in top 
quark cross section is reported as the systematic uncertainty. 

There are shape uncertainties for those spectra obtained from simulation: ti and elec- 



24 



TABLE III: Systematic uncertainties for the ti cross section are combined assuming different 
sources are uncorrelated, but shape and acceptance systematics from each individual source is 
100% correlated. 
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±3.0% 


other acceptance syst. 




±2.0% 


±2.0% 


non-W (rate) 


±1.5% 




±1.5% 


ti shape 


±1.5% 




±1.5% 


tb (single top production) 


±0.5% 




±0.5% 


total 


+12.4% 
-13.9% 


±12.3% 


+20.6% 
-21.5% 



troweak tb production. The uncertainty due to the ti shape comes from the difference 
between PYTHIA and HERWIG jj] simulations. The theoretical electroweak single top 
quark tb production cross section uncertainty is small, which is known to approximately 
3% [13]. We conservatively apply an uncertainty of 30% to the single top cross section. A 
background shape uncertainty also results from the uncertainty in the size of the ti con- 
tamination in the taggable but untagged sample and is negligible as we discussed above. 
The shape difference between the mistag and H^±HF, shown in Fig. El is small compared 
to other systematics. 
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TABLE IV: The top mass dependence of the measured total tt cross section. The acceptance and 
the leading jet Et shape depend on the top quark mass. 



mass (GeV/c 2 ) 


170 


175 


180 


cross section (pb) 


+1.6+1.3 
°-^-1.7-1.4 


6 n +1 - 5+1 - 2 

°- u -1.6-1.3 


c R +1.4+l.l 
°-°-1.5-1.2 



VII. CONCLUSION 

In conclusion, we have measured the tt production cross section in the lepton plus $t plus 
jets channel. W+ > 3 jets events were selected with at least one jet identified as a b- 
jet through secondary vertex reconstruction. Signal and background were separated us- 
ing the shape of the leading jet Et distribution. The measured total tt cross section is 
6.0^i;Q(stat.)l}J(syst.) pb where we have assumed a top quark mass of 175 GeV/c 2 jl^ . 
This is consistent with the Standard Model prediction Il4| and with the recent result from 
CDF in the dilepton channel {l-^L^.ukl ^ 0.4 pb) |2j. The measured cross section depends 
on the top mass since a heavier top produces more energetic jets. This affects both the 
signal-background shape discrimination and the acceptance. A change in the top mass of ± 
5 GeV/c 2 (12| alters the cross section by 6-8% as shown in Table ITVl 

This result also demonstrates that kinematic determination of the signal fraction using 
the leading jet Et provides good signal to background discrimination. This technique can 
be used as an effective constraint in future tt measurements, such as the top quark mass. 
This method reduces the sensitivity to statistical fluctuation of the background because the 
signal-to-noise ratio is determined from the top sample itself. 
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